Direct measurement of the biphoton Wigner function through two-photon interference 
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We show that a simple modification in the well-known Hong-Ou-Mandel (HOM) experiment 
provides the direct measurement the Wigner function of the biphoton state generated in parametric 
down conversion. In this approach, a negative value of the Wigner function is a sufficient condition 
for non-gaussian entanglement in the two photon system. In the general case, the Wigner function 
provides all the required information to infer entanglement using well known necessary and sufficient 
criteria. We analyze our results using two examples of parametric down conversion processes taken 
from recent experiments. The present work offers a new vision of the HOM experiment that further 
develops its possibilities to realize fundamental tests of quantum mechanics involving decoherence 
and entanglement using simple optical set-ups. 
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Entangled photon pairs play undoubtedly a central role 
in quantum information processing and quantum com- 
munication. Photons are the most efficient quantum in- 
formation carriers, not only for their intrinsic propaga- 
tion speed, but also for the variety of degrees of free- 
dom they possess, both discrete and continuous. Some 
examples of quantum information protocols that have 
been realized with photons are teleportation [1 , quan- 
tum key distribution [2 , one-way quantum computing 
[3] and quantum repeaters [4]. Moreover, entangled pho- 
ton pairs enable the realization of fundamental tests of 
quantum mechanics, as Bell type inequalities [5], since 
the no-signalling condition is relatively easily fulfilled. 

Spontaneous parametric down conversion (SPDC) is 
the most widely used process to generate entanglement 
in different (independent) degrees of freedom of a photon 
pair. Detecting, characterizing and manipulating this en- 
tanglement is a key issue for quantum information appli- 
cations. This problem is fundamentally different if one is 
dealing with discrete degrees of freedom (e.g. polariza- 
tion), or with continuous ones (e.g. spatial or spectral). 
While for two qubit states in the discrete case and for 
gaussian states in the continuous one, necessary and suf- 
ficient conditions exist for entanglement detection, for 
higher dimensions or more general configurations, solu- 
tions are subspace dependent [6]. However, using high 
dimensional systems and non Gaussian states leads to a 
number of important and interesting applications, such 
as entanglement distillation [7 , quantum computation [8] 
and high precision measurement |5]. For these reasons, 
understanding and classifying such states is a matter of 
importance and fundamental interest. 

Photons produced by SPDC can be highly non- 
separable because the characteristics of the pump beam 
and of the nonlinear medium are transferred to global 



degrees of freedom of the photon pair. For this reason 
we often speak of biphoton states. For instance, using 
a continuous wave (cw) pump and considering degener- 
ate, monochromatic and polarized fields, the two-photon 
state can be written as [TOl ITT] 



^+(Pi +P2)F- (pi -p 2 )|pi,P2)dpidp 2 , (1) 



where F + is the normalized momentum distribution of 
the pump beam, F_ is the phase matching function and 
Pi the Transverse Momentum (TM) vector of the i-th 
photon [10]. Eq. can indeed be obtained for sev- 
eral types of continuous variables (namely spatial or fre- 
quency coordinates) and in a wide range of experimental 
setups, as we will show below. We denote p± = pi ± P2 
and q± = qi ± q2 (the sum and differences of position 
coordinates). The correlations of the biphoton are deter- 
mined by the functions F±(p±) and their Fourier trans- 
forms .F±(q±), which describe the photons in the trans- 
verse position coordinate. To gain information about the 
entanglement in state (HI), measurement of the coinci- 
dence distributions in at least Fourier conjugate planes is 
required. Furthermore, to gain total information about 
the quantum state through tomography requires mea- 
surement of correlations along additional directions in 
the phase space of the transverse spatial variables [T2] . 
It is obvious from Eq. u\l that all of the information 
about correlations between photons 1 and 2 can be ob- 
tained through measurements on the sum and difference 
variables, p±. 

In the present paper, we show that the biphoton 
Wigner function can be measured directly using an adap- 
tation of the Hong-Ou-Mandel (HOM) interferometer 
[T3] . In the following, we detail how to detect the bipho- 
ton Wigner distribution in the whole phase space and 
consequently, obtain a full characterization of the two- 



photon quantum state, by using readily available linear 
optics elements. Furthermore, we show that HOM inter- 
ference can be used as an entanglement witness for non- 
Gaussian entanglement. Negative values of the bipho- 
ton Wigner function is a sufficient condition to prove the 
presence of non-Gaussian entanglement in the system. 

Let us first recall the principles of the HOM interfer- 
ometer, shown in Fig [I] Each photon of a pair created 
from SPDC is sent through one of the two arms of an 
interferometer. They are then recombined in a 50/50 
beam splitter (BS) and detected by detectors A and/or 
B. When the photons are indistinguishable and reach the 
BS simultaneously, they bunch and follow the same path. 
Coincidence detections in detectors in A and B are thus 
less likely, and the so-called "Hong-Ou-Mandel dip" is 
observed 13 . 
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For notational simplicity, we write the state of a con- 
tinuous variable degree of freedom of the photons as 
|\I>) = JJ F(pi,p 2 )\pi,p 2 )dpidp 2 , where pi are variables 
associated to the z-th photon (i = 1,2) that propagates 
in the i-th arm of the interferometer. We impose that all 
other degrees of freedom are separable from the consid- 
ered one, which can be guaranteed using filters to make 
local projections. In order to better illustrate the main 
idea, we will first take pi as being the TM of the z-th 
photon produced in the SPDC process, and consider only 
one spatial dimension [10 . Examples of applications to 
other continuous degrees of freedom and extension to two 
dimensions will be given below. 

In the present version of the HOM experiment, we 
suppose the interferometer is calibrated and both pho- 
tons reach the BS simultaneously. We then add a posi- 
tion translation 25 to photon labeled 2 and a momentum 
translation \i to photon labeled 1. These operations are 
currently done using linear optical elements [14]. After 
both position and momentum displacements, the bipho- 
ton state, impinging in the BS is: 



FIG. 1: Scheme of the HOM-type interferometer. Devices 
represented by boxes in each arm displace the continuous v; 
able degree of freedom in conjugate spaces, so that the Wig 
function can be measured in all the points of phase space 
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After the BS, the two-photon state is given by 

|tf B5 > = 1/2 [[d Pl dp 2 F(p u p 2 )e- 2ip * 6 (\ Pl - fj) A \P2)B ~ \Pi ~ ^)b\P2)a + \pi - ^)b\P2)b - \pi - Ii)a\P2)a). (3) 



We will now focus on the coincidence detections only, i.e., consider only states corresponding to two photons exiting 
in different paths, A and B. The coincidence probability I(n, 5) thus reads 



/ (M^)=2 _ 2 Re 



F(pi +/i,p 2 )F*(p 2 , Pl +fi)e- 2i ^-^ s d Pl dp 2 
I 



(4) 



Let us already note that separability of the biphoton 
wavefunction F(p 1 ,p 2 ) = fi(pi)f 2 (p 2 ) implies I(/jL,S) < 
1/2. This applies to two-photon mixed states as well, 
since they can be constructed as a convex sum of sepa- 
rable pure states [15]. Thus /(/i, S) < 1/2 is an entangle- 
ment witness for general two-photon state of the TM. A 
similar result was obtained for other degrees of freedom 
inRefs. [IgHIg]. 

Let us now turn to the main result of our paper. For 
the sake of simplicity, we will independently discuss each 
transverse axis, x, parallel to the incidence plane, and y, 
orthogonal to the incidence plane. This is necessary since 
there is a fundamental difference between x and y axes 
under reflection upon a vertical mirror since p x — »> — p x 
and p y — >• p y , which has interesting consequences on our 



results, as will be seen below [15] . 

We start by assuming, as in Eq. (fiT), F(pi^,p 2: i) = 
F-(p-j)F+(p+j), where i = x,y. This assumption is 
verified in most experiments with SPDC and is commonly 
used when studying entanglement in this process [10J [20] . 
From now on, in order to simplify the notation, we will 
index functions instead of variables so that F+(p+^) = 
i^+(p+), for instance. 

The integral in Eq. Q reads differently for x or y 
coordinates [T5] , For the y axis, we have 

J\F y+ (fi + p+)\ 2 dp + x (5) 

F y _(n+p_)F;_(Lt-P-)e- 2ip - 5 dp-, 



I 



while for the x axis we have 

|F x _(/i + p_)| 2 dp_ x 



(6) 



J F x+ {n + p + )F* x+ (p - p + )e- 2i ^ s dp + . 



First, let us notice that normalization can be chosen so 
that the first integrals in Eqs. ^ and ([6| are unity 
(integrals in p+ and p- , respectively) . Eq. (pi becomes 

/ 'F y _(p+p_)F;_(Lt-P-)e- 2ip - 5 dp-=TrW v -(Li,S), 

(7) 
while Eq. Q becomes: 



[f x+ (p+ P+ )F* + (p- P+ ) 



-2ip+<5 



dp + =irW x+ (p,5). 

(8) 

where W y -(ja^S) and W^+ (/i, S) are, by definition, the 
Wigner functions at point (/i, 5) [21] associated to wave 
functions F y _ or F x+ , respectively. 

Thus, the coincidence probability in this adapted ver- 
sion of the HOM experiment reveals the Wigner function 
at phase space point (/i, 5): 



I(p,5) 



;Wj(p,5), 



(9) 



where j = x+ or j = y—. It is a straightforward calcu- 
lation to show that the Wigner function of a non-pure 
state can also be directly measured using the HOM set- 
up described above [15] . 

From the Wigner function one can infer all the neces- 
sary information about the state, and in particular, en- 
tanglement for Gaussian and non-Gaussian states. The 
witness defined by i(/i, £) < 1/2 allows to detect non- 
Gaussian states since they may have negative values of 
the Wigner function. Though the witness itself does not 
detect gaussian entanglement, we can nevertheless use 
the Wigner function to test Gaussian entanglement us- 
ing other criteria [22] . These facts, added to the one 
that no assumption is being made on the width of the 
distribution, are clear advantages of the present method 
over discretization based techniques for detecting entan- 
glement in continuous variable systems [20 . 

We note that the Wigner function appearing in ([7]) is 
a single-party Wigner function referring to the sum or 
difference coordinates of the biphoton. This is a direct 
consequence of the form of state (fil) and momentum con- 
servation. 

We have shown that the reflexion asymmetry of the 
TM correlates the measurement of the entanglement 
properties of F + or F_ to orthogonal traverse directions 
(Eqs. (|7|) and (J8|). However, one can measure these 
functions in either axis, since they can be controllably 
interchanged, for instance by adding a Dove prism ori- 
entated at 45° in both arms of the HOM interferometer 
that rotates the fields by 90°. 



Up to now, we have independently considered each de- 
gree of freedom of the biphoton, but the spatial variables 
are inherently two-dimensional (2D), as exemplified by 
Eq. (fiT). This leads to a four-dimensional Wigner func- 
tion, instead of a two-dimensional one. Using the ob- 
tained results and considering the reflection properties 
of the BS, it is straightforward to show that the four- 
dimensional Wigner function returns information about 
F+ in the x direction, and about F_ in the y direction. If 
we consider, for instance, a two-photon state of the form 
(fiT), in the approximation where both transverse coordi- 
nates are separable, the application of our results to the 
two dimensions simultaneously, gives [T5] 



J- \l^xi *x'-) l^yi Vy) 



1 7T 2 



W X +(fJL X ,5 X )Wy-(lly,5y). (10) 



Using Dove prisms can lead to similar expressions as (10) 



involving orthogonal coordinates, as mentioned above. 
It is also worth mentioning that even if transverse co- 
ordinates are not separable, I(/j, x ,5 x ', /j, y ,5 y ) reveals the 



Wigner function of each transverse coordinate as in (10), 



with the difference that, in this case, we have non pure 
states in each coordinate [T5] , 

We now illustrate our results by studying in more de- 
tails some examples. Let us first consider entanglement 
in TM naturally produced in cw SPDC using a gaussian 
pump. In this case, we have, in ([5]), that F + (p + ) = 

^ - <^), where »„ 



and F- (p_ 



smc 



is the width of the pumping beam momentum distribu- 
tion, k is the wave number of the pumping beam and L 
is the non linear medium's length. For simplifying rea- 
sons, we will only study coordinate y. This can be done 
by fixing fi x and S x in x and scanning only the y phase 
space, i.e., varying fi y and 5 y only. Thus, 7rW x +(fi x ,5 x ) 
is a multiplicative constant. In the studied example, 
since the pump is gaussian, this constant is necessar- 
ily positive, and can be set to 1 by a proper choice of 
/ia, S x and the width of the pump. Thus, Eq. ( [l0| ) di- 
rectly provides W y -(fJL y ,5 y ). Supposing that spatial co- 
ordinates are separable, we can apply Eq. ^ to com- 

Corre- 



pute it, with F y -(fj, =bp_) 



J(»±p-) 2 L\ 



tsincf^ 



sponding results are shown in Fig. pa) and (|2p), showing 
Wy-(fJLy,S y ) and I(fi y ,Sy). Using realistic parameters, 
we can have I(jjl x , S x ) = 0.56 for S x ^0.1 mm, which fits 
well in the width of the transverse position distribution 
of the photon pairs. The relative violation of the separa- 
bility threshold is of over 10 %. Usually, this function is 
approximated by a Gaussian, and the relatively high vio- 
lation of the entanglement witness shows the limitations 
of this approximation. 

Our results can be further exploited by modifying 
the pumping configurations and measuring the biphoton 
Wigner function that depends on the pump profile, i.e., 
W x +(/jL x ,5 x ). We can, for instance, create Schrodinger 



cats [23] in the TM space and directly probe their Wigner 
function. This can be done by coherently splitting the 
pumping beam in two and displacing one with respect 
to the other in momentum space, for instance with the 
help of Spatial Light Modulators. As a consequence, the 
TM distribution will be centered in two different points, 
which distance we denote by Ap p . The biphoton Wigner 
function in his case is as in Fig. (J2J;) , where we considered 
Ap p = b^k/L. The entangled biphoton is highly non- 
Gaussian and violates the proposed witness by over 80% 
(see Fig. d2[i)). This configuration can also be useful to 
study the decoherence of entangled non- Gaussian states 
through the Wigner function [24] . 

Let us remark that previous results, as in [IT] can now 
be re-interpreted with the present formulation. We see 
that in [IT], as in the usual HOM experiment, the Wigner 
function at point 7(0, 0; 0, 0) was measured (see [15] for a 
revision of these results using the present formulation) . 

We provide now a second example of an experimen- 
tal set-up where our results can be applied. It consists 
in SPDC generated from a pulsed pump in semi con- 
ductor waveguides [25j [26]. In this case, pairs of pho- 
tons entangled in frequency are created, so the space 
of CV is one dimensional only. The output wave func- 



/2ttw v 



and F_(cj_) = 



tion is such that F + (co> + ) = 

J where rii are the refractive in- 



(ni-ri2)L 



smc 



f {ni—ri2)Loj- 
2V2c 011± " V 2V2c 

dex of the medium for the i-th photon, L is the medium's 
length, uu± = uo\ ± 0U2 where uoi is the i-th photon fre- 
quency. There is a complete analogy between these func- 
tions and the one dimensional TM case, and the same 
reasoning can be applied leading to the measurement of 
the Wigner function. However, scanning the whole phase 
space in this case demands using optical elements lead- 
ing to frequency displacements \i^ and frequency pro- 
portional dephasing S T . Frequency displacements can 
be realized using techniques as the one demonstrated in 
[27] , while S T displacements can be done either by time- 
delaying one arm of the interferometer or by using linear 
optics elements. Expected results with the considered 
functions are depicted in Figs. pfe), pf))- In [C]? a 
state analogous to the Schrodinger cat in Fig. (|2p), was 
created by pumping the medium with a Gaussian beam in 
a regime where two different phase matching conditions 
apply. A violation of less than 4% was observable using 
displacements in the S T axis only by time delaying one 
photon (equivalent to displacements in S x in Figs, pfc), 
(|2]i)). Negative points were observed which we can now 
interpret as interference fringes of the Wigner function 
of a biphoton Schrodinger cat state (see [15] for further 
discussion). 

In conclusion, we have provided a new interpretation 
of the coincidence probability in the HOM experiment 
in terms of a biphoton Wigner function. We adapted 
the HOM set-up so as all the points of the Wigner func- 
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FIG. 2: Wigner function (normalized to it) and region 
I(/j,i,5i) > 1/2 in some chosen variable (i = x,y, cj, r) 
for three output wave functions of the bi-photon in the 
SPDC process. Variables in all the plots are in units of 

the relevant physical parameters: fi x \JkL~ 1 , 5 X Vk~ 1 L\ , 

/io; [cL -1 (ni — 7i2) _1 ] and 5 T \c~ 1 L(n\ —712)] • (a) Output 
of a cw pumping creating entanglement in the transverse 
momentum (TM) distribution. Violation of about 12% can 
be obtained with displacements in the TM axis, as shown 
in (b). (c) and (d) Schrodinger cat state in the TM with 
Ap p = by^k/L. We see that a violation of over 0.9 is ob- 
tained, (e) and (f): Frequency entangled states produced in 
the SPDC process through pulsed pumping. 



tion can be detected by using available linear optics el- 
ements. Negative points of the Wigner function are as- 
sociated to non- Gaussian entanglement, that can be di- 
rectly detected in this type of experiment. For Gaussian 
entanglement, the information provided by the Wigner 
function can also be used applying other entanglement 
criteria. We analyze previous experimental results using 
the new perspective provided by our formulation. Our 
results open the path to realizing new fundamental tests 
of quantum mechanics in a simple and currently used 
optical setup. 
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